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Abstract 

We predict inelastic light scattering spectra from electron collective excitations in a coaxial 
quantum well embedded in a core-multishell GaAs/AlGaAs nanowire. The complex composition, 
the hexagonal cross section and the remote doping of typical samples are explicitly included, and 
the free electron gas is obtained by a DFT approach. Inelastic light scattering cross sections 
due to charge and spin collective excitations belonging to quasi-ID and quasi-2D states, which 
coexist in such radial heterostructures, are predicted in the non-resonant approximation from a 
fully three-dimensional multi-subband TDDFT formalism. We show that collective excitations can 
be classihed in azimuthal, radial and longitudinal excitations, according to the associated density 
fluctuations, and we suggest that their character can be exposed by specihc spectral dispersion of 
inelastic light scattering along different planes of the heterostructure. 
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I. INTRODUCTION 


Inelastic light scattering (ILS) spectroscopy is a fundamental optical technique to study 
electronic properties of excess carriers in semiconductorsFirst, it enables to probe elec¬ 
tronic elementary excitations, which in semiconductors occur in the range of few tens of meV, 
with near-infrared and visible light for which high-quality tunable lasers and detectors exist. 
Second, ILS enables to study separately different types of elementary excitations thanks to 
the polarization selection rules.^ Indeed, the ILS cross section arises from collective charge 
density excitations (CDEs), or plasmons, when the polarizations of the incident and scat¬ 
tered light are parallel (polarized conhguration), and from collective spin density excitations 
(SDEs) when the two polarizations are orthogonal (depolarized conhguration). In addition, 
single-particle excitations (SPEs) can be observed in both polarization conhgurations under 
strong inter-band resonance conditions.^ Finally, by adjusting the relative direction of the 
incident and scattered photons, it is possible to tune the momentum transferred to the elec¬ 
tronic system in a given direction, which, due to the momentum dependent selection rules,^ 
enables to map the dispersion of the elementary excitations. 

The wavevector dispersion of elementary excitations is strongly dependent on the di¬ 
mensionality of the electronic system.HEEl This is particularly important for intra-subband 
plasmons, as well as the Landau damping in SPEs continua. Such features have been demon¬ 
strated in ILS experiments in semiconductor heterostructures of reduced dimensionality such 
as quantum wells (QWs),®^ quantum wire^^*^^*^ and quantum dots.*^^*^ The comparison 
of the experimental results with theoretical models yielded information about the carriers 
energy structure, density, mobility and many-body contributions. 

Recently, ILS experiments and theoretical predictions have been used to conhrm the accu¬ 
mulation of excess electrons in the conduction band of radial heterostructures in nanowires, 
namely, core-multisheh nanowires (CSNW) hosting a coaxial quantum well (coQW), similar 
to the one shown in Fig. [^a).^ The extra signals observed in the ILS spectra of remotely- 
doped coQWs, in comparison with those of undoped reference samples, were shown to origi¬ 
nate from the collective excitations of a high-mobility electron gas conhned inside the coQW. 
Furthermore, ILS resonances could be selectively assigned to quasi-one-dimensional (qlD) 
and quasi-two-dimensional (q2D) states which coexist in the structure, due to the hexagonal 
nanowire cross section, being localized at the corners and facets of the coQW, respectively. 
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Such a spatial-dependent dimensionality and the tubular shape of the electron channel, 
which bridges between q2D and a qlD electron gas, guarantee a particular dispersion of the 
collective excitations in the conhned in-plane directions and a diameter-dependent momen¬ 
tum dispersion along the invariant longitudinal direction. However, results in Ref. |T3 were 
restricted to a backscattering geometry with incidence normal to the facets of the coQW 
and did not probe the energy dispersion of such excitations. 

It is important to note that ILS in coQWs substantially differs from ILS in planar struc¬ 
tures. On the one hand, in coQWs there is only one translational invariant direction (along 
the nanowire axis), the other direction being wrapped around the CSNW core, removing the 
continuous energy dispersion. The discrete symmetry also induces a inhomogeneous carrier 
distribution, with coexisting qlD and q2D states (see Fig. [^b)). On the other hand, the 
photon held always impinges on more than one coQW planes, and with different angles, at 
the same time (see Fig. [^a)). Since selection rules depend on the photon held direction, ILS 
spectra in coQWs are expected to be substantially more complex than in traditional QWs. 

The dispersion of plasmon modes has been theoretically investigated in cylindrical tubu¬ 
lar geometries.^I^^ These studies, mostly RPA calculations, illustrated the dimensional 
crossover from q2D to qlD electron gas for very small radius of the tubule, a regime appro¬ 
priate, e.g., to carbon nanotubes but not to CSNWs. Furthermore, these works neglected 
any discrete symmetry characteristic of the nanowire cross section. The lowest collective ex¬ 
citations of an electron gas in CSNWs, including the hexagonal symmetry of the electronic 
system, have been studied recently by our group,^ which enabled to rationalize the complex 
excitation spectra in terms of symmetry allowed excitations and symmetry selective Landau 
damping. However, the high electron density regime, relevant for ILS experiments, and the 
ensuing self-consistent held was not included. Likewise, ILS dispersion was not discussed. 

In this paper, we study theoretically the ILS spectra of a prototypical GaAs/Alo.aGao.rAs 
GSNW hosting a remotely doped coQW at large electron densities. h’^ l ^^l Assuming non¬ 
resonant conditions, the scattering cross section is obtained from the dynamic response 
functions of the excess electrons. Galculations of GDE and SDE spectra are performed in a 
3D real-space multi-subband DFT and TDDFT formalism which accounts for the complex 
geometry and composition of the structure, the self-consistent held of excess electrons, and 
the plasmon-phonon coupling. Three diherent classes of collective excitations arise in the 
spectra, which involve density huctuations occurring in the azimuthal, radial or longitudinal 
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directions of the coQW. We show that these collective modes can be singled out by properly 
designed ILS experiments. 


II. THEORETICAL MODEL 

A. Self-consistent ground-state calculation 

To obtain the electronic states of a CSNW we employ a standard envelope function 
approach in a single parabolic band approximation. Electron-electron interactions have 
been treated at the mean-held level through a typical Kohn-Sham (KS) LDA procedure. We 
consider a CSNW which is spatially invariant along the NW axis direction z (see Fig. Oa) for 
axis dehnitions). Therefore, we factorize the envelope functions as 

with parabolic energy dispersion Snik^) = Cn + In the r = {x^y) plane, the NW cross 
section is hexagonal and the material and doping modulations are mapped on a hexagonal 
domain using a symmetry-preserving triangular grid. The self-consistent potential vks{i^) = 
v{r) +v^{r) +v^^{r) includes the effect of the spatial conhnement v{r) determined by the 
materials conduction band offset, the Hartree potential v^{r) generated by the free electrons 
and the static doping, and an approximate exchange-correlation potential Further 

details can be found in Ref. [2S1 


B. ILS spectra calculation 

ILS spectra have been calculated within a multi-subband TDLDA non-resonant formal- 
ism.l^ The scattering cross section due to CDEs is obtained from the imaginary part of the 
momentum-dependent density-density response function 
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where, q is the in-plane and the in-wire, longitudinal components of the total momentum 
Q = (q, exchanged in the scattering process, i.e., Q = Q* — Qs, Q* and Qs being the 
momenta of the incident and scattered photons, respectively, and a; = ca* — is the ILS 
energy shift. To obtain the dynamic DDRF, we expand n(r, r', u) in terms of the in-plane 
KS envelope functions. 
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Here, Yljjqz^u) are the elements of the DDRF for the KS system, 
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2n u - {ej{kz + qz) - ei{kz)) + ig' 
where g = 2 accounts for electron spin degeneracy, fn{kz) is the Fermi occupation function 
and g is the electron damping parameter. VijkniQz) and are the direct Coulomb and 
exchange-correlation matrix elements, respectively, which describe the dynamic interaction 
of two electrons, one of which gets scattered from state i to j and the other from k to n, 
with an exchange of momentum q^. 

The direct Coulomb matrix elements read 


Vijkniqz) = J dr j dr'(j)i{r)(l)*{r)V{r-r',qz)(j)l{r')(j)n{r'), (5) 

where V{r — r', q^) is the Fourier transform of the Coulomb operator in the direction. We 
have used a frequency dependent dielectric constant entering the denominator of the 
Coulomb operator instead of the semiconductor high-frequency dielectric constant Soo, that 
phenomenologically accounts for the coupling of CDEs with the phonons of the underlying 
crystal lattice. 


-ulo + iVphU: 

Here, ujlo and ujto are the resonant frequencies of the longitudinal and transverse phonons, 
respectively, and Vph is a damping parameter. 

The exchange-correlation matrix elements read 
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with fxci'f'y ^0 being the dynamic exchange-correlation kernel dehned in the adiabatic LDA 
as the derivative of the static exchange-correlation potential with respect to the ground state 
density, fxcir,r') = - r'). 

The TDLDA formalism also provides the ILS intensity due to SDEs. These are only cou¬ 
pled through indirect-exchange interactions. Consequently, they always appear redshifted 
with respect to their plasmonic counterparts. The ILS cross section due to SDEs is therefore 
obtained from the imaginary part of the so-called momentum-dependent irreducible response 
function (IRE), n(Q,a;), which is calculated as above for the DDRF, but with the direct 
Coulomb integrals set to zero in Eq. ([^. 

From the DDRF and the IRE we can obtain the density fluctuation induced by the electro¬ 
magnetic held at a given energy and momentum, the so-called induced density distribution 
(IDD), from Kubo’s correlation formula.^ Thus, for CDEs the IDD is given by 

SncDE{r,<i,q^,uj) = j dr', ( 8 ) 

whereas for SDEs, SnsoEi'^, Q, Qz,^) is equivalently calculated substituting Il(r, r', g^,a;) by 
U{r,r',q„u;). 

III. NUMERICAL RESULTS 

A. Electron gas ground-state distribution 

Calculations are performed for the CSNW shown in Fig. [^a). A GaAs core is surrounded 
by Alo.aGao.yAs/GaAs/Alo.aGao.yAs/GaAs layers with thickness indicated in the hgure. A 
n-type h-doping layer is included in the middle of the outer Alo.aGao.rAs layer. System pa¬ 
rameters are similar to those of the sample recently studied in Ref. HTJ although the thickness 
of the GaAs coQW is slightly narrower here (26 nm in the experiment) to spectrally separate 
different classes of transitions, which are partially overlapping in Ref. [T^ as discussed in 

sec. urm 

Material parameters for GaAs (Alo. 3 Gao. 7 As) are the isotropic electron effective mass 
m* = 0.067 (0.092) and the static dielectric constant £ = 13.18 (12.24). The GaAs/Alo.sGao.rAs 
conduction band offset is taken as 284 meV. The Fermi level is pinned at the middle of the 
GaAs band gap (taken as 1.43 eV), the temperature is 4 K, and a common homogeneous 
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FIG. 1. (a) Schematics of a prototypical CSNW investigated in the paper. Dark, middle and light 
green colors are used for Alo.sGao.yAs, GaAs, and n-doping regions, respectively, (b) Calculated 
electron gas distribution over the CSNW cross section. Only the GaAs core and coQW regions are 
shown. 

triangular grid with ~ 3.21 nodes/nm^ is used to numerically integrate the KS and Poisson 
equations. Unless otherwise stated, we have used a constant density of static donors over 
the (5-doping layer of 2 x 10^®cm“^, which yields an accumulated electron gas with a total 
density ~ 1.13 x lO^cm”^. As shown in Fig. [^b), electrons are almost exclusively localized 
in the GaAs coQW with only a minority of them distributed over the GaAs core. In the 
coQW, in turn, electrons are preferentially localized in the hexagonal corners in order to 
minimize electron-electron interactions and single-particle conhnement energies. This is an 
electron distribution previously observed in different hexagonal GSNW heterostructures at 
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sufficiently large electron density,which points to coexisting channels of different 
dimensionalities, namely, qlD states in the corners and q2D states in the facets of the 
coQW. 


B. ILS at normal incidence 

We now study the ILS spectra due to the elementary excitations of the conduction free 
electrons accumulated in the CSNW in a back-scattering geometry (i.e., with antiparallel 
incident and scattered photons) and with the photon held at normal incidence with respect 
to the top facet of the coQW. In this conhguration no momentum is transferred to the 
electronic system along the invariant CSNW axis direction. In the calculation, we have 
assumed a typical excitation energy of 1.92 eV. We approximate |Qi| = |Qs|, yielding a 
transferred momentum |Q| = 1.38 x 10® cm“^, which in our sample corresponds to ~ 1.03 kp, 
kp being the Fermi wave-vector. We have used a basis set formed by the 100 lowest-lying KS 
subbands which assures convergence within the considered energy range. Finally, in order 
to account for the plasmon-phonon coupling effect we have used the GaAs high-frequency 
dielectric constant Soo = 10.86, the GaAs phonon resonances ujto(lo) = 33.72 (36.69) meV, 
and a damping parameter Tph = 1 meV. 

In Fig. we show the calculated GDE (a) and SDE (b) ILS spectra with an electron 
damping rj = 2 meV, which leads to peak bandwidths comparable to those measured in 
recent ILS experiments in GSNWs.^^ We also show a breakdown of the spectra in contribu¬ 
tions due to excitations in the GaAs coQW (blue) and core (red). As can be seen, the cross 
section is almost entirely originated in coQW excitations, with some weak core excitations 
at low energies that hardly influence the total spectrum. Likewise, we have found no appre¬ 
ciable coupling between core and coQW excitations by comparing the spectra in Fig. |^with 
equivalent ones obtained by deactivating the dynamic Goulomb and exchange-correlation 
integrals between core and coQW states (not shown). Such a result contrast with the strong 
inter-tubule coupling observed in coaxial tubuled^^^ in which, however, the inter-tubule 
separation was of the order of few A. 

Both GDE and SDE spectra are dominated by a low-energy broad band, labeled c and 
f in Fig. 1^ a) and (b), peaking at ~ 18.9 and ~ 10.8 meV, respectively. From the energy 
of the corresponding peak in the SPE spectrum (not shown), we estimate a depolarization 








FIG. 2. CDE (a) and SDE (b) spectra in a back-scattering configuration with the photons at 
normal incidence on the top CSNW facet. Blue and red lines show the contribution to the total 
spectra from excitations occurring in the GaAs coQW and core, respectively. The dashed line in (a) 
shows the GDE spectrum calculated neglecting the plasmon-phonon coupling. Vertical gray lines 
mark the position of the GaAs LO and TO phonon resonances, (c)-(g) Real (left) and imaginary 
(right) parts of the IDDs calculated at the resonances indicated by arrows in the spectra (a) and 
(b). 

shift ~ 5.8 meV and an excitonic correction ~ 2.8 meV. As reported in Ref. [T71 this broad 
signal is dne to excitations between the discrete snbbands which arise from the discretization 
of the 2D continnnm of an eqnivalent flat QW dne to azimnthal periodicity of the coQW. 
Note that in a planar QW these collective modes wonld correspond to gennine intra-snbband 
excitations and wonld not be excited at normal incidence. In onr CSNW these modes are 
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activated through the lateral facets which offer oblique orientation to the incident photons. 
By this reason, and in order to avoid confusion with the genuine ID intra-subband excitations 
along the CSNW axis, we will refer to these excitations as 2D intra-subband excitations. 

This assignment is conhrmed by the calculated IDDs shown in panels (c) and (f) of Fig. 
respectively. These clearly correspond to density fluctuations running along the lateral facets 
of the coQW with nodal surfaces normal to the coQW planes (azimuthal nodes). No nodal 
surfaces parallel to the QW plane (radial nodes) are present, since these are associated to 
inter-subband excitations lying at larger energies, as we discuss next. 

The CDE spectrum in Fig. a) also shows two less intense peaks below and above the 
GaAs phonon resonances (gray vertical lines), labeled d and e. The spectrum calculated 
neglecting the plasmon-phonon coupling (dashed line in Fig. [^a)) demonstrates that these 
two peaks are the split pair of a single CDE resonance due to coupling with the phonons. 
The corresponding SDE peak, labeled g in Fig. [^b), is at ~ 31 meV. This collective mode 
corresponds to the hrst inter-subband excitation in the coQW, as clearly demonstrated by 
the calculated IDDs shown in Figs. id) ,(e) and (g) featuring one radial node in the top 
and bottom facets. The depolarization shift for this mode (calculated using the energy 
of the CDE without phonon coupling) is ~ 3.9 meV and the excitonic correction is ~ 2 
meV. Therefore, many-electron corrections are smaller for the inter-subband than for the 
2D intra-subband mode, which is reasonable due to the higher conhnement regime of the 
former. Note that the position of the inter-subband peaks is very sensitive to the coQW 
thickness, and for thicker coQWs (as in Ref. fTTl) these modes are likely to merge to the 
broad 2D intra-subband peak. 

In Fig. [^we show the effect of reducing the value of the electron damping, which might 
be achieved in samples with higher electron mobilities, on the calculated ITS spectra. In 
contrast to a flat QW, where the damping only modihes the bandwidth of single intra- or 
inter-subband excitations, in a CSNW a small damping exposes a complex hne-structure. 
Indeed, at the energy of the broad 2D intra-subband band, numerous close by collective 
excitations show up. As discussed above, these excitations are associated to transitions 
between states with different number of azimuthal nodes, arising from the azimuthal dis¬ 
cretization of the 2D continuum. Also the inter-subband peak, which is mainly associated 
to transitions with a change in one radial node, shows several resonances at small damping, 
which correspond to different changes in the number of azimuthal nodes, as shown by the 
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FIG. 3. (a) CDE and (b) SDE spectra calculated using an electron damping of 2, 1 and 0.5 meV 
from top to bottom lines. For clarity, the curves have been offset and scaled, as indicated by labels. 

corresponding IDDs (Figs. |^d),(e) and (g)). Note indeed that the two types of excitations 
cannot be excited individually due to the different impinging angle on different facets. Since 
in Fig. [^the hne structure emerges already at a damping of ~ 1 meV, the hne structure is 
likely to be observed in future ILS experiments. 


C. Longitudinal dispersion 

We now study the dispersion of the ILS spectra as the direction of the photons is tilted 
(always assuming a back-scattering geometry) in the longitudinal direction, along the CSNW 
axis. We indicate with 6 the angle between the direction of the photon beam and the y-axis, 
and with ip the angle between the projection of the beam on the {x, z) plane with the z-axis 
(see inset in Fig. |^b)). We consider the dispersion over the {y,z) plane {p = 0°) as the 
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beam is tilted from normal to parallel to the CSNW top facet (i.e., from 6 = 0° to 6 = 90°). 
In this conhguration the momentum transferred along the invariant direction excites genuine 
intra-subband transitions within the ID parabolic subbands. At large angles, moreover, SPE 
continua gain spectral weight and induce Landau damping, which in a hexagonal NW is a 
symmetry restrictive process.^ 



ILS shift (meV) 


FIG. 4. CDE (a) and SDE (b) spectra dispersion over the {y, z) plane. Calculated spectra with 
different photon orientations 0 in steps of 10° are shown. Each curve is shifted proportionally to 
the in-wire transferred momentum q^. The top curves correspond to qz = 0, the bottom curves 
correspond to q =~ 1.03kp. The inset in (b) illustrates the scattering geometry. 

The {y, z) plane dispersion of the CDE and SDE spectra are respectively shown in Figs. 
(a) and (b). Clearly, the spectra are strongly affected by the photon orientation. A main 
feature is the disappearance of the inter-subband resonances as the photons are tilted toward 
the longitudinal direction. For SDE, e.g, the peak at ~ 31 meV disappears at 0 > 30°. The 
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FIG. 5. (a) CDE spectra dispersion over the (y, z) plane calculated using an electron damping of 
0.1 meV. Curves corresponding to different photon orientations are shown for selected values of 9 
and have been shifted proportionally to the in-wire transferred momentum q^- (b) IDD calculated 
at the resonant position of the long-lived CDE marked with an arrow in Eig. (a) for 9 = 30°. 

inter-subband phonon-coupled peaks in the CDE spectrum also disappear, but a narrow 
peak appears at ~ 39 meV as a consequence of the coupling between the phonons and the 
persisting lower energy broad band which approaches the phonon resonances as 9 increases. 

Vanishing of the inter-subband resonances with increasing transferred momentum in the 
invariant direction has been reported for planar systems both in simulation^ and experi- 
ments.^^SEZI This is due, on the one hand, to Landau damping of some collective modes inside 
SPE continua and, on the other hand, to the reduction of the momentum transferred normal 
to the QW plane. In our calculations this results from the ILS selection rules arising from 
the integrals in Eq. Q, which vanish for inter-subband transitions when q = 0 (6* = 90°). 
Correspondingly, the 2D intra-subband excitations associated with azimuthal density fluc¬ 
tuations are also expected to weaken as 9 is increased since the momentum transferred along 
the azimuthal direction is reduced. In fact, the the 2D intra-subband CDE resonance at 
~ 19 meV for 0 = 0° (top curve in Fig. |^a)) reduces in intensity as 9 increases. However, 
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an additional peak starts to be observed at 9 = 30° as a low-energy slionlder of the 2D 
intra-snbband resonance. This new peak originates in ID intra-snbband excitations and it 
shows the characteristic strong blneshift with of an intra-snbband plasmon.^^ 

In order to better analyze this resonance, in Fig. j^a) we examine its dispersion in detail 
nsing a smaller electron damping (0.1 meV). The hgnre shows that the resonance is origi¬ 
nated in a single collective mode which becomes the only one optically active for 6 > 20°. 
Therefore, contrary to the nnmerons 2D intra-snbband and inter-snbband modes showed 
in Fig. only a single intra-snbband plasmon is observed in spite of the high nnmber of 
occnpied snbbands. A similar resnlt was reported in Ref. HSI for a zero-thickness cylindrical 
tnbnle. Specihcally, when mnltiple ID snbbands are occnpied, one among all possible intra- 
snbband plasmons gets an increasingly larger spectral weight. This is the intra-snbband 
plasmon corresponding to a coherent oscillation of all electrons along the longitndinal di¬ 
rection of system. The IDDs calcnlated for the long-lived peak at 9 = 30° (Figs, [^b)) 
conhrms this assignment, since almost no density flnctnation is fonnd over the CSNW cross 
section becanse the ID intra-snbband character of the mode indnces density flnctnations in 
the longitndinal direction. 

Finally, note that in the longitndinal limit (9 > 60° in Fig. |^a)) the ILS intensity 
is exclusively due to the intra-snbband SPE continuum, and the long-lived intra-snbband 
plasmon is Landau damped. This is a consequence of the effective dimensionality of the 
system induced by the large density of states along the facets, as pure ID plasmons would 
not be Landau damped.!^ 


D. In-plane dispersion 

We next discuss the dispersion when the photon beam is tilted in the plane of the cross- 
section of the CSNW, that is the {x, y) plane (see Fig. [^b) inset). To estimate the expected 
anisotropy, we compare the spectra calculated at normal incidence (0 = 0°) and with photons 
propagating parallel to a maximal diameter {9 = 30°). ILS selection rules were discussed 
in Ref. |23] for these scattering conhgurations in a model coQW system with few occupied 
snbbands.^ While some of the excitations are excited in both conhgurations, some other 
excitations are selectively totally or strongly suppressed in either conhgurations. Similar 
arguments apply to the present study. However, in a realistic experimental situation, which 
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FIG. 6. (a) CDE and (b) SDE spectra for a sample with doping density 2 x 10^® cm“®. Black-solid 
(red-dashed) lines correspond to (/? = 90°, 9 = 0 ° {9 = 30°). The inset illustrates the scattering 
geometry, (c) CDE and (d) SDE spectra with doping density 2.7 x 10^® cm“®. (e) and (f) are IDDs 
calculated at the CDE energies marked with arrows in panel (c). 

we simulate here using a large electron gas density and electron damping in the meV range, 
it would not be possible to single out specihc excitations, due to the large density of states 
of the coQW. Nevertheless, these selection rules are at the origin of the results discussed 
below. 
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The CDE and SDE spectra for the sample studied in the previous sections are shown 
in Figs. 1^ (a) and (b). Tilting the photon beam over the {x^y) plane slightly reshapes the 
peaks but otherwise it has a small effect on the spectra. In general, at 6* = 30° the intensity 
of the 2D intra-subband peaks is increased, while inter-subband peaks are weakened, in 
comparison with the peaks at 6^ = 0°. This is originated in the fact that for 6^ = 0° more 
momentum is transfered in the radial than in the azimuthal direction, thus favoring inter¬ 
subband excitations, whereas the opposite holds for 6 = 30°. 

In Figs. I^c) and (d) we show the spectra for a larger doping density, 2.7 x 10^®cm“^, 
with a total free electron density ~ 3.11 x 10^cm“^.l2Sl Compared with the spectra for the 
lower density, all peaks appear blueshifted due to the larger subband splitting. For CDEs, 
however, this is also due to the dynamic many-electron Coulomb contributions typical of 
larger electron densities.^ Moreover, the 2D intra-subband and the inter-subband peaks 
partially overlap due to the larger blueshift of the former. In this larger carrier density 
regime, the spectra are strongly anisotropic in the (x, y) plane. For instance, in the SDE 
spectrum for 0 = 0° the inter-subband peak at ~ 31 is stronger than the 2D intra-subband 
resonance at ~ 16 meV, while the opposite is true at 0 = 30° Similar, but less pronounced, 
reshaping is also predicted for the CDE spectrum. Indeed, the IDD calculated at the energy 
of the most intense CDE peak for 0 = 0° (Fig. |^e)) shows mainly a radial density fluctuation 
in the top and bottom facets, while the most intense peak for 0 = 30° (Fig. [^f)) is rather 
originated in azimuthal density fluctuations. 


IV. SUMMARY AND CONCLUSIONS 

CSNWs are quite complex systems from the electronic point of view. On the one hand, 
due to the tubular shape, the electron gas bridges between a q2D and a qlD system. On the 
other hand, the prismatic cross-section induces an inhomogeneous distribution of carriers, 
forming qlD and q2D channels at the corners and facets of the heterostructure, respec¬ 
tively. Furthermore, in ITS experiments the photon field inevitably impinges with different 
directions on the different facets of the embedded coQW at the same time, exciting differ¬ 
ent types of excitations simultaneously. Therefore, ILS spectra in CSNWs are considerably 
more complex than in planar heterostructures. At the same time, ILS cross-section is an 
extremely informative probe of the nature of the electronic states, particularly if dispersion 
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is measured in appropriate directions, as we have discussed in this paper. 

To show the potential of ITS experiments in CSNWs, we have predicted the spectra by 
CDEs and SDEs at experimentally relevant regimes. The real-space 3D multi-subband DFT 
and TDDFT formalism used here proved of predictive quality in a previous theoretical- 
experimental study which has been conducted at normal incidence.^ Here, we have gener¬ 
alized this approach to predict the dispersion of the ILS resonances with the photon held 
rotated in different directions. We have shown that the ILS spectra hold information on 
three different type of collective excitations which form in CSNWs, namely, 2D-hke intra¬ 
subband, inter-subband and intra-subband excitations. These are associated with density 
huctuations in the azimuthal, radial and longitudinal directions, respectively, and can be 
singled out in ILS from their spectral position and bandwidth, and from their specihc dis¬ 
persion in properly designed experiments. 

In particular, i) the 2D intra-subband and inter-subband peaks vanish as the photon 
held is rotated toward the CSNW axis; ii) in this regime a single long-lived intra-subband 
plasmon can be observed; iii) the ILS spectra are anisotropic with respect to the azimuthal 
angle of the photon held, the degree of anisotropy being larger for higher electron densities. 
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